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Abstract 

We present a next-to-leading order calculation of jet production in 7*7 colli- 
sions from e + e~ scattering in a region where the virtuality Q 2 of the probing 
virtual photon is small compared to the transverse jet energy. The calcula- 
tion is based on the phase-space slicing method. The initial state singularity 
of the virtual photon is factorized into the structure function of the virtual 
photon, using the MS factorization scheme for virtual photons. Numerical 
results are presented for LEP2 conditions. The perturbative stability of the 
pure direct virtual photon approach is compared to that of including resolved 
virtual photons in different regions of Q 2 . We make predictions for cross sec- 
tions which suggest that different parametrizations of virtual photon parton 
densities should be distinguishable by measurements of jet cross sections at 
LEP. 
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1 Introduction 



Considerable progress has recently been made in investigating the structure of the virtual 
photon in jet production from eP-scattering at HERA. On the theoretical side calcula- 
tions are available in leading order (LO) and next-to- leading order (NLO) |2|, [|. An 
increasing number of experimental data becomes available || and the confrontation of 
these results with the NLO calculations indicates that the concept of a resolved virtual 
photon, i.e., the virtual photon as a source of quarks and gluons, is necessary to describe 
the data in the low Q 2 region, where Q 2 is the virtuality of the photon. The parton distri- 
bution function (PDF) of the virtual photon is constructed in analogy to that of the real 
photon, with the difference that the virtual photon PDF has an extra Q 2 - dependence built 
in, in addition to the usual factorization scale dependence. In the limiting case Q 2 — > 0, 
the virtual photon PDF's reproduce the real photon PDF's. Since so far only limited 
data exist on the structure of the virtual photon ||, the modeling of the Q 2 -behaviour 
of the virtual photon PDF is still rather ambiguous. Two LO parametrizations of the 
virtual photon PDF are available that fit the data ||, namely those of Gliick, Reya and 
Stratmann J7JI, which have very recently been updated by Gliick, Reya and Schienbein || 
(GRS), and those of Schuler and Sjostrand |J (SaS). The GRS group has also calculated 
the virtual photon PDF's in NLO, but no parametrization is available up to now, since 
the differences between the LO and NLO parametrizations are small and the available 
data is not yet very precise. 

It is desirable to find alternative ways of probing the virtual photon PDF's, either to 
test their universality or to further constrain their Q 2 - dependence. One possibility is 
the reaction 7*(Q 2 ) + 7(P 2 — 0) —>■ jets(Pr) + AT, where the virtuality of the probing 
photon Q 2 has to be sufficiently small in comparison with the transverse jet energy Ej^ 
to allow for a hadronic component in the virtual photon. The real photon with virtuality 
P 2 ~ involved in the collision has both a direct pointlike and a resolved hadronic 
part. This reaction can be obtained at e + e~ colliders by single-tag experiments, and 
preliminary results have been reported from the OPAL collaboration at LEP |10[ for 



6 < Q 2 < 30 GeV 2 . For the large Q 2 region we have recently evaluated the NLO QCD 
corrections to this process and found them to be small and to improve the scale dependence 
[ I I]. Moving towards the region of small virtualities, especially for Q 2 E^, one expects 
logarithms from the initial state 7* — > qq splitting of the form ln(Q 2 /s), where \fs is 
the partonic center-of-mass energy, to become large and spoil the convergence of the 
perturbative expansion. A procedure for subtracting these terms from the direct virtual 
photon cross sections and absorbing them into the PDF of the resolved photon has been 
worked out for the case of eP-scattering in in the framework of the phase-space-slicing 
method in analogy to the subtraction of 1/e poles from photoproduction JT^, 13fl . In 



section 2 of this paper we will work out the subtraction procedure for the 7*7-scattering 
case, closely following the approach in ||. With the subtraction performed, one has to 
include a resolved virtual photon contribution into the jet cross sections. In the deep- 
inelastic described in [fLlJI , the virtual photon couples only directly and the real 
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Figure 1: The different components contributing to jet production in 7*7 -scattering. 



photon can have a direct and a resolved component. These are called the direct (D) 
and single-resolved (SR) contributions. In addition to these, the so-called single- virtual 
resolved (SRS) and double-resolved (DR) contributions occur, when the virtual resolved 
photon components are included. The four contributions are shown in Fig. [jj In the SRS 
case, the real photon interacts directly and the virtual photon is resolved. For this case 



the NLO matrix elements can be taken from photoproduction [[E|, [13], [L4J]. Finally, in the 
DR case both photons are resolved, for which the matrix elements are described in fT3|| . 
Since NLO calculations for the case of 77-scattering with both photons being on-shell, 
i.e., Q 2 ~ P 2 ~ 0, are available in the literature [I3L O, TE\, we compare our results in 
the limiting case Q 2 — > with these calculations. 

It has been shown in |7], || [| that the parton content of virtual photons is not solely 
described by purely perturbative contributions in the region Aq CD < € /1 2 , where 
fi 2 is the probing scale, in contrast to the expectations from ||lq| . However, as has been 
pointed out in M, O, it is not clear for which values of Q 2 (and fi 2 ) the non-perturbative 



part is relevant and down to which value of Q 2 one should trust perturbation theory. In 
section 3 we will therefore numerically study the perturbative stability of the 7*7 — > jets 
cross sections in the region Aq CD <C Q 2 <,E^, where E T is the typical hard scale in jet 
production. We study the K factors for different scales and virtualities and compare the 
scale dependences of the unsubtracted jet cross sections D and SR with those obtained 
after the subtraction of the logarithms and inclusion of resolved components SRS and 
DR. From this comparison we can deduce the regions of photon virtuality where fixed 
order perturbation theory gives reliable predictions for jet cross sections and where, on 
the other hand, the parton content of a virtual photon is relevant. 

Finally, in section 4, we make predictions for low Q 2 jet production including a virtual 
resolved photon component as functions of Q 2 and E? and discuss whether the resolved 
component can be observed at LEP energies. Further numerical results can be found in 
P|. The paper ends with a summary of the results and conclusions. 
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2 Low Q 2 jet cross sections 



Taking over the conventions given in []IT] , the process we are interested in can be written 



as 



e + (k a ) + e~(k b ) ^ e + (k' a ) + e-(k' b )+jets(E T )+X , (1) 

with the subprocess J*(q a ) + l(qb) —* jets + X, where q a = k a — k' a and qb = k b — k' b . The 
virtualities are given by Q 2 = —q 2 and P 2 = —q 2 = 0. The cross section for the process 
(H) at large Q 2 is given by the convolution 



da P 



dQ 2 dy a dy b 



J2 J dx b F j/e -(y b )f b/j (x b ) 
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where y a = (q a k b ) / \k a k b ) . The variable y b G [0,1] describes the momentum fraction of 
the real photon in the electron. The momentum fraction of a parton in the real photon 
is x a G [0, 1] and the PDF of the real photon is f b /-y(x b ). The indices U and L denote 
the unpolarized and longitudinally polarized virtual photon contributions. Finally, the 



function F 1 / e - (y b ) is the Weizsacker- Williams approximation for the real photon [TJ5| 
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where m e is the electron mass and 9 max is the maximum scattering angle of the untagged 
electron. An improved version of the equivalent-photon approximation in e + e~ collisions 

leading however only to small corrections. We use the simpler 



m 



has been derived in 

formula since our studies are exploratory and we do not compare with data. 

For the SR case, the 
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The definition of the partonic cross sections da^l b are given in 
NLO calculations can be taken from [2T|, whereas the NLO calculations for the D case 
have been presented in |llj; both calculations employ the phase-space-slicing method to 
extract the singular phase-space regions of the real corrections. The singular integrals of 
both the real and the virtual corrections are handled in dimensional regularization. Most 
of the l/e n poles for the real and virtual corrections cancel and the remaining initial state 
corrections are subtracted into the real photon PDF. We are in the following interested 



in the initial state splitting for the virtual photon, 7* 



QQ: 



which has been evaluated in 



for the SR contributions. Here, we repeat this calculation for the D process^ 

i*{q a ) + i(qb) -> g(pi) + q(P2) + g(ps) ■ 



(4) 



The first step is to extract from the 2 — ^ 3 matrix elements of the reaction ([|) a term 
with the characteristic denominator from the 7* — > qq splitting which gives the singu- 
lar contribution in the limit Q 2 — > 0. We call this term Hk- In the same singular 

1 Note, that in Q the photon virtuality corresponding to our Q 2 was called P 2 . 
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limit the three-body phase space of the qqg final state, dPS 1 ^, factorizes according to 
dPS (3) = dPS (2) cffS (r) , with dPS (2) being the usual two-body phase-space and dPS (r) 
being the phase-space of the the singular region [fj, [13], ^TJ. We define the variable 
z a = (P2P3) I (q a qb) £ [r]a, 1] that gives the fraction of the momentum q a that partici- 
pates in the subprocess after a particle has been radiated in the initial state. The variable 
r] a G [0, 1] is connected to z a through r\ a = x a z a . The term H K is integrated over the 
singular phase space up to a cut-off y s with the result 

[dPS^H K = ^(ArraQ 2 ) 2 f^2N c C F M(Q 2 )T 1 (s,t,u) , (5) 
J 2,71 J z a 



where 



M(Q 2 ) = —P q ^{ Za ) In 1 + (6) 



is the singular term and 



Xy{s,t,u) = -- - - (7) 
u s 



is the LO photon-parton scattering matrix element where s, t and u are the usual Mandel- 
stam variables. The photon splitting function is given by P q ^(z a ) = Nq (z 2 + (1 — z a ) 2 ). 
The term @ is large for Q 2 C s and singular for Q 2 = 0, as expected. It is the same 
universal term as obtained in for the eP-case. We therefore introduce, in accordance 
with 0, the subtraction term 

/ M 2 \ 

r q ^(z a , M 2 ) = In [ Q2{1 "_ Za) J P^r(za) - N c (8) 

which is to be absorbed into the PDF of the virtual photon. After this subtraction the 
remaining finite term in M(Q 2 ) yields 

M (Q 2 )-m - w >* ( j^nj ) +5 • w 

In addition to the singular term ln(M^ / Q 2 ) two finite terms have been subtracted in order 
to achieve the MS factorization for Q 2 7^ [[|. It is defined by the requirement that the 
remaining finite term M(Q 2 )jj^ is equal to the finite term obtained after factorization of 



the real photon initial state singularities, which can be found in [T2|, |M[ . Of course, 
this has consequences concerning the selection of the PDF of the virtual photon, the de- 
tails of which can be found in [0. This completes the calculation of the contribution from 
the virtual photon initial state singularity. Note that the above described subtraction 
only concerns the transversely polarized virtual photons in reaction (|]), since the con- 
tributions for longitudinal photons vanish in the limit Q 2 — > 0. Furthermore, up to now 
no virtual photon structure function for longitudinal photons has been constructed. The 
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final formula for jet cross sections including resolved virtual photon contributions in e + e 
scattering at low Q 2 is a generalization of (p]) and reads 



^g^- = Zjdx a dx b F l/e -(y b )f bh (x b )^ 



Va 



-rj r (x a ,Q 2 )da ab 



, 2(1 y a ) . \jL 
H 5(1 - x a )do ■ , b 

Va 



(10) 



where fa l fy*( x a,Q 2 ) is the virtual photon PDF and x a G [0,1] describes the momentum 
fraction of the parton in the virtual photon. The direct photon interactions are included 
in this formula through delta functions. For the direct virtual photon one has the relation 
f™ r /^*da ab = 5(1 — x a )da^, b , whereas for the direct real photon the relation is / 7 / 7 <icr a & = 
5(1 — x b )da lb) where da ab refers to the partonic cross section. 

Formula ( |TUD is implemented in the fixed higher order program JetViP [^2|| , with which 
all numerical results in this paper are produced. The input parameters for all numerical 
studies are the following. We assume LEP2 conditions, i.e., the energies of the incoming 
leptons are E a = E b = 91.5 GeV. We integrate over the full range of y a ,y b G [0, 1] and 
use the value 9 max = 0.025 in @. For the real photon we always employ the PDF's of 
Gliick, Reya and Vogt (GRV) P5fl , whereas for the virtual photon we will use the new 
GRS || parametrization and the SaSID PDF || transformed to the MS-scheme @. The 
GRS PDF's are constructed for Np = 3 flavours, the production of the heavier c and 
b quarks is supposed to be added as predicted by fixed order perturbation theory. We 
calculate our cross sections with Np = 5 flavours and use the two-loop formula for the 

(5) 

strong coupling constant without threshold effects with Aq<^ d = 153 MeV, even for the 
LO results. The cross sections are plotted in the 7*7 cms, where also jets are defined. 



We use the Snowmass accord [24], where two partons i and j are recombined, if for both 



partons the condition R it j < R is fulfilled, where R it j = y (Vi ~ Vj) 2 + (4>i ~ 0j) 2 an d 
rjj, (fij are the rapidity and the azimuthal angle of the combined jet respectively, defined 
as 

E Tj = E Tl +E T2 (11) 

E Tj T]j = E Tl 7]i + Et 2 T] 2 , (12) 

E Tj 4>j = E Tl 4>i + E T2 4> 2 ■ (13) 

We choose R — 1. 

As a numerical check of our calculations and of the consistency of the MS-scheme for 
virtual photons, we compare in Fig. ||| a,b the cross section (|T0|) in the limit Q 2 — > with 
existing calculations for photon-photon scattering WM. On the virtual photon side, the 



Q 2 is integrated out by using the Weizsacker- Williams approximation [1{J with Q 2 max = 1 
GeV 2 . In Fig. |2] a the comparison is made for da/dEp as a function of Et, where the 
rapidity has been integrated out in the range 77 G [—2,2]. The dots are the results from 



Kleinwort and Kramer fl4|l , whereas the curves are the predictions from JetViP, for all 
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four components, D, SR, SRS and DR, as discussed above, now using the SaS PDF's. 
One sees a perfect agreement. The SR and SRS distributions are very close to each 
other, but do not coincide exactly, since different PDF's are employed for the real and the 
virtual photon. However, the SRS and SR cross sections agree exactly on the partonic 
level. The excellent agreement seen in the Et distributions holds also for the rapidity 
distribution da/dr] of Fig. b for all four components. The transverse energy in these 
curves has been integrated out with Et > 3 GeV. The SR and SRS components have their 
respective maxima at opposite sides of the rapidity range. The SR distribution is peaked 
at positive rfs, since the direction of the virtual photon was taken to be the positive z-axis. 
In photoproduction, the DR component gives considerable contributions in the low Et 
region, as can be seen from the Et spectra. The D cross section only starts to dominate 
for Et larger than 10 GeV. The SR and SRS components are of minor importance. 



3 Direct vs. resolved virtual photon approaches 

We now compare the two approaches for calculating jet cross sections at small Q 2 , namely 
considering the direct coupling of the virtual photon only and, on the other hand, including 
the resolved virtual photon components. We cover the region Aq CD C Q 2 ^ E t where 
one could expect contributions from resolved virtual photons to be of importance. For 
the discussion in this section we will only use one of the virtual photon PDF's, namely 
those of SaS, which are built for Np = 5 flavours. The input parameters for the following 
plots are as described in the previous section. 

In Fig. ^ we have plotted the ratio of the single inclusive jet cross sections 

where da NLO contains all Born terms. The cross sections are obtained by integrating 
over the rapidity-range |?7| < 2 and for four different values of transverse energy, Et — 
3, 7, 10 and 25 GeV, in the region Q 2 e [0.1, 200] GeV 2 , where the lower limit is roughly 
Q 2 min — 10 x Aq CD . The renormalization and factorization scales [ir and \ip are set equal 
to Et- We observe that K is close to unity for all four curves for Q 2 ^ 30 GeV 2 which 
indicates the perturbative stability of the cross sections in the deep-inelastic region. The 
ratio exceeds 1 below Q 2 = E\ and rises monotonically towards smaller Q 2 for the three 
smaller E^s. For Et = 25 GeV the NLO terms give only very small corrections below 
5% in the whole Q 2 region down to Q 2 min - For E T = 3 GeV the rise towards smaller Q 2 is 
strongest and leads to a NLO correction of nearly 100 % near Q 2 ~ Qmin> but is already 
around 50 % for Q 2 ~ E T /5. For the larger Et values the NLO corrections do not rise 
that dramatically, however below Q 2 = 10 GeV 2 the corrections exceed 20 % for E T = 7 
and 10 GeV. 

We take out the case E T = 3 GeV to make a more detailed study of the perturbative 
stability by looking at the scale dependence of the cross sections. In Fig. || a-d we show 
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the single jet inclusive cross section integrated over \r)\ < 2 and Et > 3 GeV as a function 
of fi/E T £ [|,3], where \i = [Ir — //p, i-e., renormalization and factorization scales 
are varied at the same time. The curves are shown in the four bins of virtuality Q 2 £ 
[0.25,0.5], [1,2], [2,5] and [5, 10] GeV 2 . Since we are interested only in comparing direct 
and resolved virtual photon contributions, we have introduced the following combinations 
of the four components D, SR, SRS and DR: the virtual direct (VDIR) is the sum of 
D and SR components, whereas the virtual resolved (VRES) is the sum of SRS and 
DR components. The component labeled VDIR^ is the virtual direct component after 
subtraction of the term (§) in the D and the corresponding term in the SR component, 
to be found in 0. 

The scale variation of the NLO VDIR contribution between the smallest and largest /x 
value in Fig. £| a amounts 30% and goes down to around 15% in the largest Q 2 bin. 
We have checked that the renormalization scale dependence alone gives a variation of 
60% for the smallest and 25% for the largest Q 2 bin, i.e., the scale variation behaviour 
becomes stable and well behaved for Q 2 approaching E T . Most of the variation with the 
renormalization scale stems from the SR component, whereas the D component is much 
more stable and varies at most 10%. This is understandable, since the D component 
is only 0(a s ) in the strong coupling, whereas the SR component is 0(a 2 ). The rather 
strong variation of the VDIR component at small Q 2 values together with the large K 
factor observed above, indicates the need for a resummed approach at small Q 2 . We have 
therefore also plotted the NLO VRES component in Fig. f| a-d. To avoid double counting, 
the VRES component can only be added after the perturbative 7* — > qq splitting has been 
subtracted from VDIR (giving VDIR5) since these terms are contained in the pointlike 
part of the virtual photon PDF. One sees two effects of including the VRES component. 
First, the scale variation is considerably reduced in the first two Q 2 bins, namely by a 
factor of 2 in the first bin and even a factor of 4 in the second bin. For the third bin the 
scale variation of the SUM = VDIR5 + VRES is more or less the same as for the VDIR 
component alone, i.e., around 20%. For Q 2 £ [5, 10] GeV 2 the SUM prediction becomes 
unstable for h<^0.7Et, since in this region Q 2 > fi 2 and the contribution from the virtual 
photon PDF is very small. But even for \i > Et the scale variation of SUM is much larger 
than for the VDIR alone. The second observation is that including the VRES component 
in NLO gives a relatively large correction to the pure VDIR component of 20-30%, up to 
40% in bin c. 

To disentangle the effects of resummation and non-perturbative parts in the virtual photon 
PDF from those of higher order contributions in the matrix elements we plot in Fig. [5] 
a-d the same curves as in Fig. [| a-d with the VRES component in LO only. The sum 
of NLO VDIR5 and LO VRES gives rather small corrections to the pure VDIR of about 
15% in the smallest Q 2 -bin and only 5% in bin b. In the bins c and d the two approaches, 
NLO VDIR and NLO VDIR5+LO VRES, give nearly the same results. The conclusion 
from this is that the virtual photon PDF is to a very large extent given by the splitting 
term for Q 2 > 2 GeV 2 at Et — 3 GeV, i.e., for Q 2 ^ |/i 2 . At smaller virtualities one sees 
effects from resummation and some non-perturbative input. Only small improvements 
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can be seen from including the VRES component in the bins a and b with respect to the 
scale variation, in contrast to the findings of Fig. |j. This is however not surprising, since 
the VRES matrix elements are only LO. Including the VRES parts in NLO reduces the 
scale dependence, as we have seen above, for bins a and b, i.e., for Q 2 ^ |/x 2 . We have also 
looked at scale variations of cross sections with larger minimum transverse energy. We find 
the overall improvement of including a resolved virtual photon component increasingly 
smaller, the larger the P^'s are. The upper bound Q 2 ^\[i 2 for the improvement of the 
scale dependence is also found for these larger P^'s. As an important result, we found no 
improvement for Q 2 > 10 GeV 2 , no matter how large the was. This is also supported 
by the Q 2 -fndependent small K factor for E? = 25 GeV in Fig. |3|. These findings confirm 
the expectations from Jl6j that for large enough photon virtualities one should end up 



with a fully perturbative prediction irrespective of the probing scale. 

The bottom line from these observations is that the resolved virtual photon approach 
improves the scale dependence und thus the perturbative stability of the jet cross sections 
for Q 2 < 10 GeV 2 as long as Q 2 <> \E\. It is interesting to see that this result, based on 
perturbative calculations in NLO QCD and without making use of the GRS PDF's, agrees 
with the restrictions implemented by the GRS group into their old parametrizations of the 
virtual photon parton densities [^]. Also for the new PDF's ||, for which no restrictions 
have been implemented, the authors stress that the resolved photon approach is only 
meaningful for Q 2 <C /i 2 , typically Q 2 <,j^fi 2 . Our results should have some relevance 
also for jet production at HERA, since the SR and DR contributions, which occur in 
eP-scattering, are included in this analysis and give considerable contributions. We have 
checked that the limitations we give above on the Q 2 -range for the improvement of the 
perturbative stability also hold for the SR and DR contributions alone. However, the 
perturbative stability in the eP-scattering case should be checked in more detail under 
HERA conditions. 

A point one has to keep in mind in this discussion is that the resolved, especially the 
DR, matrix elements may give important contributions to jet cross sections, even though 
the non-perturbative input from the virtual photon PDF is small. As we have seen in 
Fig. |], the sum of NLO VDIR5 and NLO VRES gives rather large corrections to the NLO 
VDIR result, even for Q 2 , where the photon PDF is given mainly by the pointlike terms. 
These NLO VRES corrections convoluted with the leading logarithmic contribution to the 
photon to quark splitting is a leading logarithmic approximation to the full NNLO result 
with a pure direct virtual photon contribution, which is so far not available. The effect 
of higher order terms in the resolved matrix elements can be much pronounced in specific 
phase-space regions, as has been shown in f25| . There, the NLO DR matrix elements are 



needed to explain the forward jet production cross section at low XBj in eP-scattering at 
HERA in the region Q 2 e [5, 100] GeV 2 . 
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4 Jet cross sections at LEP2 



We now turn to predictions of cross sections under conditions to be met at LEP2. We 
keep the numerical input as specified at the end of section 2. We start with absolute 
predictions for single jet inclusive cross sections integrated over rapidity \rj\ < 2 



as functions of the transverse energy. In Fig. |6] a-d the spectra are shown for Q 2 = ^,1,2 
and 5 GeV 2 . We have plotted all four components to see their relative importance. To 
be able to include the resolved contributions we have subtracted the logarithmic terms 
in the direct components, leading to the Ds* (dashed) and SRs (dash-dotted) curves. 
One sees the strong fall-off in Et and the decrease of the absolute values with increasing 
Q 2 . For Q 2 = 0.1 GeV 2 the DR component is rather important for the whole shown 
Ex region and especially dominates at Et = 3 GeV. This is in accordance with the 
photoproduction results presented in Fig. || a. Movin towards larger virtualities, the 
influence of the DR cross section becomes smaller, and at Q 2 = 5 GeV 2 the D component 
is the dominant contribution. Only for the smallest Et values is the DR cross section in 
the same magnitude as the D one. The SR and SRS contributions are always small. In 
addition, for the SRS cross section a further suppression can be seen for increasing Q 2 due 
to the suppression of the virtual photon PDF. The virtual resolved components have a 
stronger fall-off with Et than the virtual direct contributions and the relative importance 
of the VRES cross sections diminish with increasing Q 2 . 

To see in more detail how large the fraction of resolved contributions in the jet cross 
sections are, we have calculated the ratio of the VRES over the VDIR one-jet inclusive 
cross sections in LO integrated over \rj\ < 2 for three different values Et = 3, 5 and 7 GeV 
for different virtual photon PDF's as functions of Q 2 in the region Q 2 G [0.1,5] GeV 2 . 
Note, that in the LO case no subtraction of the logarithmic terms has to be performed. 
As a scale we have chosen fiR = = Et- The results are shown in Fig. [?| a-d. For all 
three PDF's in Fig. |7| a-c one observes the expected fall off with rising Q 2 . The fall off is 
stronger for the smaller scales. For the larger scales the VRES component diminishes with 
respect to the VDIR component, but reaches out farther into the larger Q 2 region. For 
the SaS parametrizations the ratio is around 0.8 below Q 2 = 0.5 GeV 2 for Et = 3 GeV 
and falls off to 0.5 for Et = 7 GeV. At Q 2 = 5 GeV the VRES cross section gives about 
15-25% of the VDIR one. These results hold also for the SaS2D case, shown in Fig. |7| b, 
although the SaSID and SaS2D parametrizations have rather different Q 2 -behaviour due 
to the different evolution starting scales of Qq = 0.6 GeV and Q$ = 2 GeV, respectively. 
The SaS2D decreases stronger and is flatter for larger Q 2 than the SaSID PDF. This can 
also be seen in Fig. [7| d for the scale Et = 5 GeV, where the two SaS parametrizations 
can be directly compared as the dashed and the dotted lines. In the same figure also 
the GRS parametrization || (full line) is shown, which produces a still different fall-off 
behaviour. For the smaller Q 2 the decrease is similar to the SaS2D case but continues 




(15) 
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stronger for the larger Q 2 . This behaviour is seen for the GRS PDF in Fig. [7] c for all 
three scales. For comparison we have plotted in Fig |7| d also the old version of the GRS 
PDF (dash-dotted line), which is rather similar to the new parametrization. It is clear 
from these curves, especially Fig. [7] d, that the Q 2 dependence is rather different for the 
GRS and SaS parametrizations which reflects the ambiguities due to the limited data on 
the virtual photon structure function. 

We could repeat the calculation of these ratios for the NLO case, but the VDIR and 
VRES components alone depend rather strongly on /ip. Only their sum is independent 
from ix p- It is preferable to adopt the strategy used experimentally to distinguish direct 
and resolved cross sections for dijet events. The OPAL collaboration have presented 
distributions in x^ for 77-scattering p6fl , where 



± _ Hjets(E±PzJ . 

t ~~ V (E + n v { ] 

Z-^hadrons\ J - / - 1 - f 



The sum in the numerator runs over the two largest jets in the event. The direct 
dominated cross sections then corresponds to those contributions, where x^ > 0.8 and 
the resolved dominated correspond to x^ < 0.8. The comparisons with NLO calculations 
show good agreement ||13|| . We however do not show curves like this here since it is not clear 
whether the statistics in the 7*7-case will be high enough to extract such distributions. 

Instead, we compare absolute single-jet cross sections like for the pure VDIR and 
the VDIR5+RES approaches with different virtual photon PDF's GRS and SaS, which 
are shown in Fig. [S] a-d. For Q 2 = 0.1 GeV 2 the VRES approach yields about a factor 
of 2 larger cross sections. This is expected, since this small virtuality lies in the pho- 
toproduction domain. The predictions from the SaS and the GRS PDF's give similar 
results. However, for Q 2 = 1 GeV 2 these two parametrizations already give rather dif- 
ferent predictions. The GRS curve is for the smaller Et's around 50% larger than the 
VDIR, whereas the SaS is 70% larger. Furthermore, the SaS curves do not fall off that 
strong with increasing E?- In Fig. § c, the difference between the two VRES predictions 
is quite pronounced. At Q 2 = 5 GeV 2 , the SaS curve gives still a 30-40% larger cross 
section than the VDIR, whereas the GRS curve yields basically the same result as the 
VDIR. In general, the difference between the VRES and the VDIR cross sections in the 
low Q 2 region is large enough that it should be possible to distinguish between these two 
approaches experimentally. Furthermore, the difference between the VRES curves with 
the SaS and GRS parametrizations is rather large and it should therefore also be possible 
to distinguish between these specific PDF's and to constrain the Q 2 -dependence of the 
virtual photon PDF. 



5 Summary 

We have presented a calculation of jet cross sections in 7*7 scattering from e + e~-collisions 
at low Q 2 in NLO QCD, employing the phase-space slicing method to extract singularities 
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in the real corrections. Logarithmic contributions from the virtual photon initial state 
have been subtracted and absorbed into the resolved virtual photon structure function. 
Comparison with existing photoproduction calculations in the limit Q 2 — > showed very 
good agreement. 

We have studied the perturbative stability of two approaches to low Q 2 jet production, 
namely the pure direct coupling of virtual photons and, secondly, the inclusion of re- 
solved virtual photons by looking at K factors and scale dependences. We found the 
resolved virtual photon approach to improve the perturbative stability for Q 2 < 10 GeV 2 
with Q 2 <; \E\. However, the NLO corrections to the resolved matrix elements may be 
important also for larger Q 2 . 

We further made predictions for inclusive jet production at LEP2 and found that it should 
be suitable to experimentally distinguish between different parametrizations of the virtual 
photon structure functions, since the DR contributions are rather important at low Q 2 . It 
should also be possible to restrict the Q 2 dependence of the parametrizations by comparing 
the data to the predictions for Et spectra of the jet cross sections. 
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Figure 2: Comparison of single inclusive jet cross sections in the limit Q 2 — > (lines) with 
the calculations of Kleinwort and Kramer (dots). The full lines give the DR, the dashed 
the D, the dotted the SRS and the dash-dotted the SR contributions, (a) Et distribution 
with 1 77 1 < 2; (b) rj distribution with Et > 3 GeV. 




Figure 3: Ratio of NLO to LO single jet cross sections integrated over \q\ < 2 for E T = 
3, 7, 10 and 25 GeV as a function of Q 2 . 
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Figure 4: Scale dependences of single inclusive jet cross sections integrated over < 2 
and E T >3 GeV as a function of fi/E T for (a) Q 2 G g, \\ GeV 2 , (b) Q 2 e [1,2] Ge^ ; 
^ g 2 G [2,5] GeV 2 and (e?j Q 2 G [5, 10] GeV 2 . The dashed line is the NLO VDIR which 
has to be compared with the full line, giving the sum of the subtracted virtual direct (NLO 
VDIRs, dotted line) and the NLO virtual resolved (NLO VRES, dash-dotted line). 
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Figure 5: Scale dependences of single inclusive jet cross sections with the same conditions 
as in Fig. $ a-d, only here the virtual resolved contribution is included in LO instead of 
NLO. 
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Figure 6: NLO single jet inclusive cross sections integrated over \rj\ < 2 as functions of E T 
for different virtualities , employing the SaS PDF's. (a) Q 2 = 0.1 GeV 2 ; (b) Q 2 = 1 GeV 2 ; 
(c) Q 2 = 2 GeV 2 ; (d) Q 2 = 5 GeV 2 . 
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Figure 7: LO ratio of VRES over VDIR of one- jet inclusive cross sections, integrated over 
\r)\ < 2 for three different scales Et = 3, 5 and 7 GeV for different virtual photon PDF's. 
(a) SaSID; (b) SaS2D; (c) GRS; (d) comparing the PDF's from a-c at Et = 5 GeV 
together with the old GRS PDF. 
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Figure 8: NLO single jet inclusive cross sections integrated over \q\ < 2 as functions of 
Et for Q 2 values as in Fig. a-d; the full line gives the unsubtracted direct (DIS) and 
has to be compared with the curves denoted GRS and SaS. These are the sum of NLO 
subtracted direct and NLO resolved cross sections, employing the GRS (dashed) or the 
SaS (dotted) parton densities. 
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